Quantum Nondemolition Measurement of a Kicked Qubit 
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We propose a quantum nondemolition measurement using a kicked two-state system (qubit). By 
tuning the waiting time between kicks to be the qubit oscillation period, the kicking apparatus 
performs a nondemolition measurement. While dephasing is unavoidable, the nondemolition mea- 
surement can (1) slow relaxation of diagonal density matrix elements, (2) avoid detector back-action, 
and (3) allow for a large signal-to-noise ratio. Deviations from the ideal behavior are studied by al- 
lowing for detuning of the waiting time, as well as finite-time, noisy pulses. The scheme is illustrated 
with a double-dot qubit measured by a gate-pulsed quantum point contact. 
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I. INTRODUCTION 

Motivated by general interest in mesoscopic quantum 
phenomena, as well as by specific interest in quantum 
computation, there has been intensive experimental and 
theoretical effort to study coherence properties and mea- 
surement possibilities of solid state qubits. In particular, 
there has been great interest in different quantum mea- 
surement proposals. This paper is concerned with a type 
of weak measurement, referred to as a Quantum Nonde- 
molition (QND) measurement. QND measurements are 
important when one wishes to extract information about 
a quantum state over a period of time without destroying 
it. This type of measurement has been proposed and im- 
plemented in optics, ^'^ gravity- wave antennas,'^ and vari- 
ous realizations of harmonic oscillators*. QND measure- 
ments have also been proposed as a quantum computing 
measurement method for single spinsi^ 

A weak measurement of a quantum system occurs 
when the state of the system is not perfectly correlated 
with the detector, in contrast to the usual projective 
measurement. Weak measurement has an advantage over 
projective measurement in that one can study the coher- 
ent oscillations of a single quantum state over a long time 
period, whereas with projective measurements, one must 
repeat the measurement many times to see the quantum 
coherent oscillations. However, weak measurements have 
the disadvantage that they usually act back on the sys- 
tem, so that one is never sure whether the result stems 
from the original quantum state alone. QND schemes try 
to overcome this limitation by coupling the detector to 
the system via an operator that represents an approxi- 
mate constant of motion of the quantum system. In this 
way, the detector disturbs the system in a minimal way. 

Recently, Averin proposed a QND scheme for a two- 
state system that involves a time dependent coupling 
which follows the qubit in spin-space while weakly mea- 
suring it.^ The scheme was illustrated with a Josephson- 
j miction system as an example of a mesoscopic qubit. 
The purpose of this paper is to generalize Averin's qubit 
nondemolition measurement. The method we propose is 
quite general and could be implemented in many different 
detector realizations. The basic idea is that rather than 




FIG. 1: The quantum nondemolition measurement performed 
with a quantum point contact. A double-dot with charge Qi 
is capacitively coupled to a quantum point contact. Current 
through the quantum point contact is periodically switched 
on and off with the help of a gate voltage Vg with the waiting 
time between pulses given by the Rabi oscillation period of 
the double dot. 



try to chase the spin around in spin space, we let it come 
to us. This is done by strobing the qubit/detector cou- 
pling constant in time, so the detector only has access to 
the qubit at approximately discrete moments. In this re- 
duced problem, by fixing the waiting time between pulses 
to be a multiple of the qubit oscillation time Tq , the ob- 
servable we want to measure is static in time. Thus, by 
measuring over many pulses, a large (back-action free) 
signal-to-noise ratio may be obtained. The liability of 
this proposal is that one cannot detect what is happening 
while the coupling is turned off, so the qubit must have 
a reasonably high quality factor. However, this scheme 
has several advantages over the original proposal: (1) the 



coupling to the qubit needs only one qubit operator (the 
one we want to measure), (2) an arbitrary qubit Hamil- 
tonian is allowed, (3) in the ideal case, no assumptions 
need to be made about the density matrix of the qubit, 
and (4) in the non-ideal case, the pulses may be noisy. 
These advantages make this kicked QND measurement 
feasible with existing technology. 

Although kicking the qubit in time sounds rather vio- 
lent, it is in fact coherent because we choose the pulses 
to be fast, so they may be described by a series of delta 
function kicks. Dephasing is induced by the fluctuating 
measurement apparatus. The inspiration for the present 
paper is from a venerable model of quantum chaos, the 
kicked rotor. *'^ Here, the situation is even simpler be- 
cause the kicked quantum system has only two states. 
Nevertheless, it is the play-off between the waiting time 
between kicks and the system dynamics that yields the 
interesting physics. In this way, the kicked QND mea- 
surement is similar to kicked atoms which give rise to 
accelerator modesiS or quantum resonancesii However, 
the idea of "stroboscopic measurement" is actually quite 
old, and the basic idea appears already for harmonic os- 
cillators in Refs. [121111^ and has recently been proposed 
for the nanomechanical resonator^ 

An implementation of the quantum nondcmolition ex- 
periment is shown in Fig. 1. A double quantum-dot rep- 
resenting a charge qubit is capacitively coupled to a quan- 
tum point contact (QPC). The QPC is opened and closed 
with the help of a plunger gate with voltage Vq ■ Before 
considering this specific physical system, we first present 
general theoretical considerations. After the detector as- 
sumptions are reviewed in Sec. ^I ttic Hamiltonian is 
introduced and the ideal quantum nondcmolition mea- 
surement is described in Sec. IHII Deviations are then 
allowed in Sees. IIVI and Ivl with detuning between the de- 
tector pulses and the qubit oscillation time, as well as 
noisy pulses of finite width. We then describe the de- 
tails of the quantum double-dot measured by the QPC 
in Sec. El (see Fig. 1). Sec. lVHI contains our conclusions, 
and Appendix IXI contains our perturbation method. 



II. DETECTOR ASSUMPTIONS 

A realistic, effective detector has two necessary char- 
acteristics: the measurement output should behave clas- 
sically and the response time of the detector should be 
much shorter than the period of the qubit oscillations. 
As mentioned in the introduction, there can be strong 
detectors whose states are well correlated with the state 
of the system and thus make a projective measurement 
of the qubit's state, or there can be weak measurement 
which only extracts partial information about the quan- 
tum state. For the weak measurement, we assume that 
linear response applies i^ Quantitatively, if the bare de- 
tector variable is /, (which we treat as a classical random 
variable), coupled to the qubit with the operator az, and 
the detector output is O, then it is given by a noisy, 



random part q, plus a small component that reflects the 
dynamics of the qubit)iS*i^ 



Oit)=qit) + ^az{t), 



(1) 



where A is the response coefficient which indicates how 
sensitive the detector is, and az is the qubit operator that 
we are trying to measure. With these assumptions, the 
temporal correlation function of the detector's output is 



C'it) = (0(t + t)0{T))r -= 2TTSg Sit - 

K{t) = {l/2)TT[paz{t)azm, 



0) 



A^ 



K{t), 



(2) 



where Sq is the spectral weight of the noise q, the density 
matrix of the qubit is p, and we have made use of the 
fact that the detector response time is much faster than 
the qubit oscillation time. Cross terms of the form {q az) 
vanish because the detector noise is uncorrelated with the 
qubit oscillations. The spectral density of the detector is 
now given by 
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(3) 
This is the object that is experimentally accessible for 
low frequencies. If the detector is coupled to the qubit 
statically, over several oscillations the measurement de- 
vice will measure the desired quantum operator iaz ) , but 
there will be detector back-action altering the qubit dy- 
namics. This back-action imposes a fundamental limit 
on the signal-to-noise ratio of 4. — We turn now to the 
pulsed scheme to overcome this limitation. 

Instead of following the spin, we let it complete the 
Rabi oscillation and couple only after a full rotation. The 
time scale of the on/off transition should be slow to the 
detector (so linear response still applies), but fast to the 
qubit. If the effective qubit has higher energy states, then 
the time scale of the kick must also be sufficiently slow so 
as not to excite these states. The dynamic coupling will 
result in a ladder of pulses at the measurement output 
device, typically of height A. The measurement of the 
quantum state of the qubit is extracted by examining 
how the rungs of the ladder are correlated with each other 
above the white noise background. 



III. THE IDEAL SCENARIO 

We first consider the simplest case where the duration 
of the pulse is zero. The qubit is then coupled to the 
detector by periodic delta functions in time of waiting 
time r, with / describing the random detector variable. 
The Hamiltonian for this ideal scenario is given by 



rr e ^ 

H^-az + -a, 



f: Sit-nr)^+H,, (4) 
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where Hf is the detector Hamiltonian. The evolution op- 
erator corresponding to this time-dependent Hamiltonian 
may be found from 



W* = T exp 



dt'H{t') 



(5) 



where T is the time ordering operator. The time integra- 
tion gives the free propagator when the delta function is 
zero, while during the kick, we can integrate the delta 
function over a small time interval e to find the qubit 
propagator for the kick: 



lim exp < — i 



r+e 



dt' 



Hg + 5{t - T 



<yzl{t) 



= exp[-icr2/(T)/2]. 



(6) 



The contribution of the free propagator is proportional 
to e in this interval so it vanishes when we take £ — > 0, 
while the kick contribution stays finite. This allows us to 
factorize the propagator as 

U^ = UjtUk:, Ujr = exp[-^r(e/2)a^ - it{1^/2)(j^I 

n 

Z^K - expHa,/(T,)/2], U'\i:) ^\{U,\i:). (7) 

Here, n = t/r is the discrete time index. It is important 
to note that Ui changes during each time step because 
/ has its own dynamics while the coupling is turned off 
that is governed by Hf. So far, the waiting time interval 
T may be arbitrary. In the ideal case, we choose r to be 
Tq = A-K l^/e^ + A^, the oscillation period of the qubit. 
Then, the free propagator is simply the identity matrix, 
so the total propagator is composed oiUicS only. 

Now that we have the propagator, all dynamical prop- 
erties of the qubit may be investigated. The first impor- 
tant property is that of the coherence of the qubit. De- 
coherence effects are caused by the randomly changing / 
variable. If we begin kicking with an initially pure state 
t/i — (a, (3) in the z basis, the density matrix changes in 
time as 



p(n) 



Z*a*l3 



Zaf3' 

l/Ci|2 



(8) 



where Z = JliLi ^^ ~^ exp[m(/)— n {5p) /2] by the cen- 
tral limit theorem. Here, we see explicitly that the pulsed 
coupling keeps the diagonal matrix elements constant in 
time, but there is exponential dephasing controlled by 
the fiuctuations in / whose dephasing rate is simply the 
noise of the detector variable. It is clear that informa- 
tion about the x and y variables is being destroyed by the 
measuring apparatus, while we have carefully preserved 
the z information of the state. 

One simple but important quantity to calculate is the 
expectation value of Oz (n) , 



tti 



For the ideal case, the propagator is just a function of 
CTz, so az is a constant of motion. Tr[pcr2] — \a\^ — |/3p, 
so the expectation value is a measure of the population 
that is time independent. We may now calculate the time 
correlation function of the detector output. The discrete 
version of Eq. |2Jl is 



C{n) = 27r5, 5nfl + {\^ /2) Tr[pa,f7,(n)], 



(10) 



where az{n) = lA^^azhi^ ■ Again, ct^ commutes with Z-/, 
a'l — 1 and Trp = 1, so Kin) never decays with n, be- 
cause every rung on the signal ladder is identical. If we 
calculate the discrete version of the spectral function, 
Eq. (PJ, the sum diverges at oj = 0. This tells us that 
the accumulated signal about Gz may be arbitrarily larger 
than the noise of the detector, S*,, regardless of how small 
the response coefficient A is. 



IV. DETUNING 



We now generalize this ideal behavior to more re- 
alistic situations. In this section, we let the waiting 
time T be close to the qubit oscillation time, but not 
equal. This gives us a perturbation parameter 7, where 
r/Tq = 1-^7. We take 7 <C 1 and let p = (1/2)1 for cal- 
culational simplicity. Noting the identity exp[— i(a/2)n • 
1 cos a/2 



a\ = icosa/^ — ifi ■ cr sin a/2 with a = rVe^ + A^ 
and the unit vector fi chosen to give the free propaga- 
tor (jiz = Te/a, Ux — rA/a), we choose a = 47r + 7 
implying hijr = exp[— 1(7/2)72 • a]. This is the same as 
free evolution by an angle 7 of our choosing. Time evo- 
lution is now different from the ideal case because the 
free and kick evolution operators do not commute. We 
may investigate the time correlator of Gz by expanding 
hi to second order in 7 to find 



K{n) = 1 



^)' f--g-« 



(11) 



{<Jz{n))=Tr[pU'<JzU'']. 



(9) 



where the sum describes the interaction with the detector 
and the Xj are recursively defined by: 

Xi = e'f\ X, =e'-^^(l + Xj_i), j = 2,...,n-l. (12) 



The fi are chosen from the same probability distribution 
and are assumed to be independent, so we can average 
over the /'s using {e^^) = exp[i(/)-S'/2] where S = {5f^) 
is the dephasing rate found from Eq. ||SJ|. The sum in 
Eq. Hll|) may now be carried out to yield 



{K{n)) = 1 - 



A 



Tql 



2tt 



n e-(^/2)("+i)+e-^/2(n-l)-7ie-^ 
2 ^ (1 - e-S/2)2 



(13) 



With the assumption 7 <C 1, Eq. ^l'A\ contains two (di- 
mensionless) rates, the dephasing rate, S, and a much 
smaher "nondemohtion" rate. It is this nondemohtion 
rate that controls the decay for times longer than the 
initial transients, {K{n)) = 1 — Q^n, 



Gd = 



A 



Tql 



27r 



1 



=S/2 _ I 



(14) 



Although this represents only the second order contri- 
bution in 7, we can extrapolate to longer time by expo- 
nentiating, {K{n)) fa exp[—Qdn]. Including the noise of 
the detector Sq, the correlation function determines the 
spectral density S{u}) = (l/27r) ^J^__^ C(n) coswn of 
the detector output. Besides the back-action, there are 
always other relaxation processes that are not well con- 
trolled. We therefore introduce Ge as a phenomenological 
environmental relaxation rate to find 



Sito) =Sg + 



A2 



smh{Ge 



47r sin^(w/2) + sinh^(ge/2 + Gd/^) 



(15) 



The sin^(ti'/2) reflects a fundamental periodicity in the 
experiment: the effect is still there if we wait some mul- 
tiple of the oscillation period. If we further assume that 
S ^ 1 (the detector noise is weak), than we can expand 
the decay rate Gd and divide by Tq to obtain the real time 
decay rate. 



rrf = 



A7 

"2^ 



2Tq 

{5p 



(16) 



This special case is analogous to Averin's line-shape.- 
The whole point of the scenario is that Eq. (|15|l describes 
the natural line-width of the qubit with detector back- 
action as a tunable parameter. 



V. NOISY PULSES 

In experiments, the ladder of pulses will have finite 
width and their own noise associated with them. In this 
section, we take these effects into account. To model this 
realistic situation, we consider the following Hamiltonian: 



2 Z tg 



Hj. 



(17) 



The variable U{t) generalizes the delta function coupling 
f{t)6{t — tit) in Sec. IIIII by broadening it into a pulse 
that lasts for a short time tq with a large prefactor 



r = TqJTQ. While the pulse is on, the detector variable 
may be decomposed into an average height with fluctu- 
ations: C/(i) — U\- SU{t). Like (/), U simply serves to 
renormalize the a^ prefactor, e — e + rU, and will not 
enter into the final results. The fluctuations are treated 
in the white noise approximation, 

{5U{t)) = 0, {6U{t)5U{t')) = {5U^)6{t - t'). (18) 

Note that although U(t) has units of energy, the conven- 
tion of Eq. (|18|l defines {SU'^) to also have units of energy. 
Now, we wish to calculate the ensemble averaged correla- 
tor K{t), Eq. ^, for one pulse after making the diagonal 
assumption for the density matrix to simplify the analy- 
sis. The method we use is analogous to the discrete case. 
The equations of motion are straightforward. 






- + r SUit)]ay, &y = [e + r SUit)]ax - Acr^ 



Act,, 



(19) 



The time correlator Eq. ^ may be calculated by formally 
integrating Eqs. 119|l in time to a given order n, as shown 
in the appendix. The closed algebra of the Pauli matrices 
allows the solution to be formally written in terms of n 
time ordered integrals. The condition r ^ 1 selects a 
certain class of terms in the perturbation theory that 
may be averaged and summed to all orders yielding 



(Kit)) 



AH 



A" 



{SU^ 



{{SU^)r^y 



-{SU^)rH 



(20) 

The basic relaxation mechanism is the same for the de- 
tuning as it is for the noisy pulses of finite width, so 
Eq. H20|l bears a close similarity to Eq. I|13|l . Notice also 
that e is not present in the formula (to this order in per- 
turbation theory). When the pulse is turned off to let the 
qubit evolve on its own, it will be the bare energy e that 
determines the qubit time Tq. Thus, the average variable 
U does not appear in the results. There are two rates 
in Eq. ^. The first, T^ = r^ {6U'^), is an additional 
contribution to the fast dephasing rate. The second is 






(21) 



which clearly corresponds to the slow nondemohtion rate 
of the noisy pulses. The parameter r~^ = Tc/Tq plays the 
role of the detuning parameter 7 in the previous section. 
We now wish to combine all the phenomena discussed. 
When both the detuning 7 and the width of the pulses 
r~^ are small, a double nondemohtion measurement is 
taking place. In this case, we can write (in the vicinity 



of w = 0) an approximation to the spectral noise of the 
detector as a Lorentzian, 



5M = Sg + 



X' 






r = re + rd + r„, 



' A2 



((5[/2 



(22) 



Here, the total broadening width of the spectral noise, F, 
has been decomposed into contributions from the envi- 
ronment, the detuning of the pulses from resonance, and 
the finite width of the pulses. 



VI. IMPLEMENTATION WITH THE 
QUANTUM POINT CONTACT/DOUBLE DOT 

To see how this nondemolition measurement will work 
in practice, we illustrate it by using a double quan- 
tum dot (DD) as the qubit and the current through 
the Quantum Point Contact (QPC) as the measure- 
ment apparatus] ^^i^7ii?ii?i29i2ii22i2;^ Additionally, a gate 

is placed above the QPC so that one may bias the barrier 
of the electrons in the QPC (Fig. 1). In this fashion, we 
can pulse the QPC so as to let current flow through dur- 
ing a short time when the gate voltage is raised above the 
barrier of the QPC. This pulsing is similar to the rectan- 
gular voltage pulses applied in the experiment of Ref . |2J 
and theoretically described in Ref. |23 The magnitude 
of this current will be sensitive to whether the electron 
is on the upper or lower part of the DD. At the detector 
output, there will be a current pulse train that may be 
analyzed statistically. The results of the nondemolition 
measurement will be seen in the spectral density of the 
current shot noise, Eq. H22|l . 

The charge response time2& trc of a QPC is very fast. 
Here we assume that even the single particle dynamics 
of the QPC governed by tqpc = h/{eV) (where V is 
the applied bias) is fast compared to the pulse duration 
tg which is in turn fast compared to the Rabi oscilla- 
tion period of the qubit (we now restore the h for clar- 
ity). Therefore, a good separation of time scales may be 
reached for the pulse time tq, 



Tqpc '^tq 't:Tq 



(23) 



For solid state qubits, there is another constraint that 
true spins do not have: more states at higher energy. 
To avoid exciting the qubit into these states, a separate 
constraint is that tq > h/AE, where AE is the energy 
difference to the next excited state. In this limit, we can 
use the physical setup described by Pilgram and one of 
the authors,^^ The interaction Hamiltonian that couples 
the DD qubit to the detector is He = QiQ/Ci where Q 
is the charge in the vicinity of the QPC, playing the role 
of / and Qi is charge on the DD. The current operator 
to be measured takes the form of a bare current lb play- 
ing the role of q, plus a self-consistent screening current 



proportional to the the potential U of the QPC, 



i = h- xu, 



X-2^E^^M 



(24) 



The (non-equilibrium) linear-response function x is de- 
fined in terms of the (large) bias V across the QPC, and 
the transmission probabilities of channel n, Tn- To ex- 
press this in the form of Eq. |^, U must now be written 
in terms of the charge Qi on the DD and the capacitances 
defined in Fig. I, U = [Cp/{C,e'^D)]Qi, where D is the 
density of states at the Fermi energy in the scattering 
region, C~^ = C^^ + C^""^ + D~'^ is the polarization ca- 
pacitance, and Qi — euz- The deviation of the average 
current is given by 



(A/(t)) = -(a,(0), 
while the shot noise spectrum is given by 

Sii{uj) = Si,i, + —k{uj), 



(25) 



(26) 



where K{u}) is the Fourier transform of K{t). Thus, the 
response coefficient A is simply the difference in current 
output depending on whether the electron is in the upper 
or lower dot. In terms of the system parameters, it is 
given by 



A 



eCo 



dT, 



nhQD 



T-im^^n 



(27) 



Physically, A is proportional to the difference in trans- 
mission caused by the location of the DD electron, 5T = 
(dT/dE) 5E, where 6E is the inverse density of states 
D''^ and Gp/Gi is the (small) coupling constant. 

By doing an experiment where one varies the time be- 
tween pulses and plots either the width of the line-shape 
of the spectral density or the signal-to-noise ratio as a 
function of the waiting time, the point of minimum width 
(or maximum signal-to-noise) finds both the oscillation 
time Tq of the qubit as well as the system's natural (un- 
modified) relaxation rate. This is because of the fact 
that the width is simply a sum of the measurement back- 
action and uncontrolled environmental relaxation. 



VII. CONCLUSIONS 

We have proposed a novel quantum nondemolition 
measurement on two-state systems with a kicking mea- 
surement device. By tuning the waiting time between 
kicks to be the oscillation period of the qubit dynamics, 
the weak measurement process does not act back on the 
desired operator being measured, though it does destroy 
the complimentary information in the conjugate vari- 
ables. Deviations from this ideal scenario are taken into 
account by allowing for a detuning of the waiting time as 



well as for noisy pulses of finite width. These effects may 
be described by an approximate Lorentzian line-shape, 
whose width is the sum of three terms; the uncontrolled 
environmental relaxation rate, a QND detuning contri- 
bution, and a QND noisy pulse-width contribution. This 
novel QND proposal was illustrated with the quantum 
double-dot, measured by a gate-pulsed quantum point 
contact. The results of the nondemolition experiment 
may be obtained from the shot noise. 
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APPENDIX A: PERTURBATION METHOD 



We now show how to obtain Eq. l(^ from Eq. |[TH|) 
and Eq. (|19|l . considering the fact that the duration of 
the pulse is much shorter than one coherent oscillation. 
We use time-dependent perturbation theory by writing 
nested time ordered integrals to a given order, and then 
take averages over all the (5[/'s. The primary object of 
interest is the correlator K{t) = (l/2)Tr[az(i)az(0)]. At 
the initial time, Tr[cr:r(0)cr^(0)] = TT[ay{0)az{0)] = 0, so 
only the terms that involve ^^(O) will contribute to the 
correlator K{t). We thus generate integral expressions 
for az{t) that end in ^^(O). As a simple example, the 
correlator to second order in time is found from 



J 



a^) = / dt'Aay{t')+a^{0)= / dt'A 
Jo Jo 



dt" [{e + r6U{t")) a,{t") - Ac7,(t")] + a,(0) + a,(0). (Al) 



To compute K{t), we set all the time arguments of the ct's to 0, and take the trace with (Jz{Q). Thus K{t) 
1 — A^ Jq dt' J^ dt" = 1 — A^t^/2. To go to higher orders, the trick is to ehminate a^- 



ayit)=ay{0)+ I dt'\~Aa,{t')~ I dt" [{e + rSUit')) {e + rdU{t")) + A^] ay{t") \ . (A2) 

I 



The order of the expansion is given by the number of 
different time integrals. In general, there are many dif- 
ferent terms to compute at each order. However, because 
the time t is much shorter than the oscillation time, the 
dominate contribution is the one given by the maximal 
number of r's, or SU's. At order 2n, such terms are of 
the form A^r^"~^. To average over the fluctuating po- 
tentials, we use Wick's theorem to average the 2n — 2 
(5f/'s in terms of m = n — 1 pairwise contractions. We 
then compute average time ordered integrals of the form 



it's] 



{6U{ti)6U{t2)...SU{t2m)). 



(A3) 



In general, contractions of this form have every possi- 
ble time-ordering of two dU^s. However, the contrac- 



tion {SU{ti)SU{tj)) comes with a delta function enforcing 
ti — tj. If time-ordered integrals occur where ti and tj 
are not in adjacent time intervals, then this term is zero 
under the integration. The only term left will be the one 
where all time arguments are time ordered. Therefore, 
the integral Eq. jSHl is given by {6U'^)"'t"'+^ /{m + 2)!. 
We may now sum these dominate contributions at every 
order in time to obtain. 



A 2.2 oo ( ^\n 



n+l 



1- 



AH 



n=2 

A2 



(5L/2)r2 ((5L/2)r2)2 
which is the same as Eq. H2()(l . 



(^^-{SU-)rH _ ^^ ^ (A4) 
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